Direct evidence for predominantly phonon-mediated pairing in high-temperature 

superconductors 
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The spectra of the second derivative of tunneling current d 2 //d\^ 2 in the high-temperature super- 
conductors YBa2Cu307_<s and Bi2Sr2CaCu20s+5 show clear dip and peak features due to strong 
coupling to the bosonic modes mediating electron pairing. The energy positions of nearly all the 
peaks in — d 2 //dV 2 -like spectra match precisely with those in the phonon density of states obtained 
by inelastic neutron scattering. The results demonstrate that the bosonic modes mediating the 
electron pairing are phonons and that high-temperature superconductivity should arise primarily 
from strong coupling to multiple phonon modes. 



The identification of phonon anomalies in the tunnel- 
ing spectra [l[ and the observation of an isotope effect [2| 
provide important and crucial clues to the understand- 
ing of the microscopic pairing mechanism of supercon- 
ductivity Q in conventional metals. For high-transition- 
temperature (high-T c ) copper-oxide superconductors, ex- 
tensive studies of various unconventional oxygen-isotope 
effects have clearly shown strong electron-phonon inter- 
actions and the existence of polarons 4^ @, H, 0, Hol . 

. Neutron scattering |l6l|. angle-resolved 
photoemission (ARPES) [l7j . and Raman scattering [l8| 
experiments also provide evidence for the strong electron- 
phonon coupling. However, these experiments do not 
provide direct evidence that pairing is predominantly 
mediated by phonons. Recent ARPES data concerning 
the effects of electron-boson interactions on electron self- 
energies have led to contradictory conclusions. On one 
hand, in the diagonal direction of momentum space (jr, 
7r), the observed 'kink' features around 70 meV in the 
band dispersion of several cuprate compounds appear to 
provide evidence for strong coupling between electrons 
and the 70 meV Cu-0 half-breathing mode [19] . The cou- 
pling to this 70 meV phonon mode is shown to be primar- 
ily responsible for d-wave high-T c superconductivity [20l ] . 
On the other hand, Dcvereaux et al. propose that the 36 
meV Big Cu-0 buckling mode rather than the 70 meV 
mode is the main player for d-wave high-T c superconduc- 
tivity based on ARPES data near the antinodal direction 
21 1 . To further complicate matters, a third group [22| 
has even provided a completely different interpretation 
for the ARPES data in the antinodal direction, that is, 
the kink feature is related to a magnetic resonance mode 
rather than a phonon mode. It is clear that the identity 
for the bosonic modes mediating the electron pairing is 
inconclusive due to these contradictory interpretations. 

In order to unambiguously identify the origin of the 
bosonic modes mediating the electron pairing, it is es- 
sential to precisely determine the energies of the bosonic 
modes. If the energies of the bosonic modes have a 
one-to-one correspondence to the energies of the peak 
positions in the phonon density of states determined 
from inelastic neutron scattering, one can definitively 
conclude that phonons predominantly mediate the pair- 



ing. For conventional superconductors, strong electron- 
phonon coupling features clearly show up in single- 
particle tunneling spectra P, (23|. The energies of the 
phonon modes coupled to electrons can be precisely de- 
termined from tunneling spectra. More specifically, the 
energy positions of the peaks in — d 2 //dV 2 , measured 
from the isotropic s-wave superconducting gap A, cor- 
respond to those of the peaks in the electron-phonon 
spectral function a 2 (uj)F(Lu) (Refs. Q, Sj|). Therefore, 
if the coupling strength a 2 (uj) does not have pronounced 
structures, the structures in the phonon density of states 
F(lo) will have a one-to-one correspondence to those in 
— d 2 //dy 2 . Such a conventional approach to identify 
the electron-boson coupling features would have been ex- 
tensively applied to high-T c superconductors if the su- 
perconducting gap were isotropic. Since the supercon- 
ducting gap is highly anisotropic in high-T c supercon- 
ductors, it is difficult to reliably determine the energies 
of bosonic modes if tunneling current is not directional. 
This may explain why the electron-boson coupling struc- 
tures extracted from earlier tunneling spectra are not re- 
producible a mong different groups [24j . |25| . |26| . In a very 
recent article 27J attempting to show an important role 
of phonons in the electron pairing, the authors assign the 
energy (52 meV) of a peak position in +d 2 I/dV 2 spectra 
of Bi2Sr2CaCu20g+5 to the energy of a phonon mode. 
Such an assignment is incorrect because the energies of 
phonon modes are equal to the energies of dip positions 
rather than peak positions in +d 2 I/dV 2 (Refs. [ll. I23]]). 

Here we report d 2 //dy 2 -like spectra for the 
high-Tc superconductors YBa2Cu307-i (YBCO) and 
Bi 2 Sr 2 CaCu208+5 (BSCCO). We find that the energy 
positions of nearly all the peaks in d 2 J/dV r2 -like spectra 
match precisely with those of the peaks in the phonon 
density of states obtained by inelastic neutron scatter- 
ing. Such excellent agreement between tunneling and 
neutron data can be explained only if the tunneling cur- 
rents in these junctions are highly directional and the 
phonon modes primarily contribute to the electron pair- 
ing. 

Figure la shows the normalized scanning tunnel- 
ing microscopic (STM) conductance spectrum g^M — 
dI/dV)/{dI/dV) B at 4.2 K for a slightly overdoped 
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FIG. 1: a) Normalized STM conductance qnm = 
(dJ /dV) / '(dl / 'AV)b versus energy E = eV for a slightly 
overdoped YBa2Cu 3 7 _ a (YBCO) with T c = 90 K, where 
(dI/dV)B is the background conductance at high biases. The 
STM conductance spectrum was taken on the (001) crystal 
face with a Pt-Ir tip at 4.2 K [28(]. b) The renormalized con- 
ductance g = gNM /Nbcs versus uj = eV — A (the energy 
measured from the gap A) for YBa2Cu307_^. The solid line 
is a smoothed curve obtained using a cubic-spline interpola- 
tion. 



YBCO with T c = 90 K, where (dI/dV) B is the back- 
ground conductance at high biases. The data are taken 
from Ref. [28| and the normalized conductance is en- 
larged by a factor of 1.06 to ensure unity at high biases. 
The spectrum shows negligible zero-bias conductance and 
sharp peak at eV =19.6 meV, suggesting that the super- 
conducting gap A is about 19.6 meV. The solid red line 
is the BCS density of states: N BC s = eV/^/(eV) 2 - A 2 , 
which lies in the middle of the data points. In Fig. lb, 
we plot renormalized conductance g — gNM /Nbcs ver- 
sus u) = eV — A (the energy measured from the gap). 
After the renormalization, the structures due to strong 
electron-boson coupling can be clearly seen in the first 
derivative spectrum of tunneling current. 

In order to precisely determine the energies of the 
bosonic modes coupled to electrons, it is essential to take 
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FIG. 2: a) — A(dg/duj) spectrum for the YBa2Cu307_a crys- 
tal together with the phonon density of states F(ui) obtained 
from inelastic neutron scattering [2S|. The left scale is for 
— A(dg/duj) and the right scale is for F(u). The vertical 
dashed lines mark peak/shoulder features in — A(dg/du). b) 
— 2A(dg /duj) spectrum together with the phonon density of 
states F(co) for the conventional phonon-mediated supercon- 
ductor Pb. Here g = N s /(N n NBCs), N s is the superconduct- 
ing density of states, and N n is the density of states in the 
normal state. The left scale is for — 2A(dg/du) and the right 
scale is for F(ui). The data of N s /N n are from Ref. PJ, and 
the data of F(lo) are from [3(| ■ The vertical dashed lines mark 
peak/shoulder features in — A(dg/dui). 



the derivative of the renormalized conductance g. In 
Fig. 2a, we show the — A(dg/dw) spectrum of the YBCO 
crystal together with the phonon density of states F(ui) 
obtained by inelastic neutron scattering 29]. For com- 
parison, in Fig. 2b we show the —A(dg/du>) spectrum 
and F(u)) for the conventional phonon-mediated super- 
conductor Pb. It is apparent that the —A(dg/dui) spec- 
tra for both YBCO and Pb show peak, dip, or shoulder 
features. A single broad peak or shoulder in F(ui) and 
— A(dg/du>) can occur when the separations of multiple 
peak features are too small compared to their widths. For 
Pb, nearly all the peaks/shoulder features in — A(dg/duj) 
match precisely with those in F(uj), as expected from 
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the phonon-mediated pairing mechanism. However, the 
peaks at 3.3 meV and 8.9 meV in the phonon density 
of states correspond to the dips in the tunneling spec- 
trum. This mismatch is due to the fact that the coupling 
strengths for the two modes are much weaker. 

If strong coupling to phonon modes also happens 
in high-T c superconductors such as YBCO, then the 
peak/shoulder features in — A(d<?/dw) should line up 
with the peak/shoulder features in the phonon density 
of states as well. Indeed, nearly all 13 peak/shoulder 
features in — A(dg/duj) of YBCO match precisely with 
those in F(u>). The strong coupling features at 7.1 meV 
and 90.8 meV in — A(dg/dui) cannot compare with these 
neutron data since the energy positions of these features 
are outside the energy range of the neutron data. Nev- 
ertheless, the phonon peak at about 6.8 meV is clearly 
seen in the high-resolution neutron data of BSCCO 31 1 
(also see Fig. 3a). The strong coupling feature at 83.7 
meV does not appear to have a correspondence in the 
neutron data where there is a strong localized incoherent 
peak at 87 meV [32^ . However, a clear peak feature at 
about 83.5 meV is seen in the earlier neutron data where 
the localized incoherent peak at 87 meV is not present 
[29l ] . The feature at 90.8 meV should be the composite 
phonon energy of 7.1 meV and 83.7 meV since the sum 
of 7.1 meV and 83.7 meV is equal to 90.8 meV. This 
is expected from the conventional strong-coupling theory 
[l|, I23I ] . Such excellent agreement between neutron and 
tunneling data provides clear evidence that the bosonic 
modes mediating the electron pairing are phonons and 
that the tunneling current of this junction is highly di- 
rectional. 

It is interesting that the low-energy modes at 7.1 
meV, 20.3 meV and 25.5 meV have much stronger cou- 
pling than the higher energy modes and thus contribute 
more to the electron pairing. In contrast, the electron- 
phonon coupling strengths for the 14.6 meV and 31.9 
meV phonon modes are much weaker because the dip 
features in — A(dg/duj) occur at these phonon energies. 

If the strong electron-phonon coupling features in 
YBCO are intrinsic, they should also show up in other 
90 K superconductors. In Fig. 3a, we show the 
— A(dg/doj) spectrum for a BSCCO crystal together 
with the phonon density of states F(ui) at 296 K. Here 
g = N s (ui)/[N n (u})NA(u})} and Na{oj) is a smoothed 
curve of N s (uj)/N n ((jj), which is proportional to 1 + 
0.99exp(-w/18.7) + 0.13 exp[-(w - 71.3) 2 /779] and does 
not have fine structures. The data of N s (w)/N n (ui) are 
obtained from Fig. 5a of Ref. [II]. The N s (uj)/N n {uj) 
spectrum has a peak at eV = 25.2 meV (Ref. .26]), sug- 
gesting A ~ 25.2 meV. This gap size is smaller than the 
average gap expected for a BSCCO with T c — 93 K. Thus, 
the spectrum may be probing a higher doping region with 
T c ~ 89 K in the intrinsically inhomogencous BSCCO. 

From Fig. 3a, we can clearly see that nearly all 
the peak/shoulder features in the — A(dg/duj) spectrum 
match precisely with those in the room-temperature 
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FIG. 3: a) — A(dg/dto) spectrum for a Bi2Sr2CaCu208+i 
crystal with T c = 93 K and the phonon density of states 
F(uj) at 296 K. The high-resolution neutron data of F(u>) 
are reproduced from Fig. la of Ref. [3lj . The left scale is 
for — 4A(dg/du) and the right scale is for F(ui). The vertical 
dashed lines mark peak/shoulder features in — 4A(dg/dui). 
b) The — d 2 (ReE)/do; 2 spectrum (E is electron self-energy) 
[33 ] and the —dg/dcu tunneling spectrum of BSCCO crys- 
tal. The — d 2 (ReE)/dc<j 2 spectrum is reproduced from 
Ref. " 
the - 



Vertical dashed lines mark the peak positions in 



-d 2 (ReE)/da; 2 spectrum. 



phonon density of states. The peak feature at 65.0 meV 
in — A(dg/du) is slightly off from the peak at 64.0 meV 
in the room-temperature neutron data. It is interesting 
to note that two peaks at 14.4 meV and 31.8 meV in 
the phonon density of states correspond to the dips in 
the tunneling spectrum. For YBCO, this mismatch oc- 
curs at 14.6 meV and 31.9 meV (see Fig. 2a), which are 
very close to those (14.4 meV and 31.8 meV) for BSCCO. 
Such quantitative agreement indicates that the 14 meV 
and 32 meV phonon modes in both YBCO and BSCCO 
have much weaker coupling to electrons. 

In Fig. 3b, we compare the strong coupling features 
revealed in the tunneling spectrum and in the electron 
self-energy spectrum of BSCCO. The spectrum of the 
second derivative of the real part of electron self-energy 
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— d 2 (Re£)/dcj 2 is reproduced from Ref. [33(. It is strik- 
ing that the peak features in — d 2 (ReE)/da-i 2 match very 
well with those in — A(dg/du>). This excellent agreement 
between the tunneling and self-energy spectra further 
demonstrates that the observed strong-coupling features 
in both spectra are intrinsic. 

It is known that any bosonic modes that mediate 
the electron pairing show up in tunneling spectra. The 
present tunneling spectra of both YBCO and BSCCO 
clearly demonstrate that the bosonic modes below 100 
meV, which contribute to the electron pairing, are only 
the phonon modes. On the other hand, pairing mecha- 
nisms based on strong coupling to a magnetic resonance 
mode predict a pronounced dip feature at about 40 meV 
in the dl/dco spectra of 90 K superconductors [13]. How- 
ever, the dip feature occurs at about 26 meV, which is 
nearly independent of doping j2?| • This suggests that 
the dip feature is not caused by strong coupling to the 
magnetic mode. Moreover, the observed large reduced 
gap (2A/kBT c > 5) implies a large coupling constant 
(>2) and a quite low logarithmic mean boson energy 
Tilj\ u (<40 meV) (23| . This rules out pairing mechanisms 
based on strong coupling to high-energy magnetic and/or 
charge excitations. Therefore, these tunneling spectra in- 
dicate that high-temperature superconductivity should 
arise primarily from strong coupling to multiple phonon 
modes. 

In summary, we have shown that the energy positions 
of nearly all the peaks in — d 2 //dV 2 -like spectra of 
YBa 2 Cu307_(5 and Bi 2 Sr2CaCu208+(5 match precisely 
with those of the peaks in the phonon density of states 
obtained by inelastic neutron scattering. The results 
demonstrate that the bosonic modes mediating the 
electron pairing are phonons and that high-temperature 
superconductivity should arise primarily from strong 
coupling to multiple phonon modes. The identification 
of an extended s-wave gap symmetry in the bulk [35| 
supports the phonon-mediated pairing mechanism. 
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